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ABSTRACT 

The pair-weighted relative velocity dispersion of galaxies provides a measure of the thermal energy 
of fluctuations of the observed galaxy distribution, but the measure is difficult to interpret and is very 
sensitive to the existence of rare, rich clusters of galaxies. Several alternative statistical procedures have 
recently been suggested to relieve these problems. We apply a variant of the object-weighted statistical 
method of Davis, Miller, & White (1997) to the Las Campanas Rcdshift Survey (LCRS), which is the 
largest and deepest existing redshift survey that is nearly fully sampled. The derived one-dimensional 
dispersion on scales ~ Mpc is quite low: ui = 126 ± 10 km s~^, with a modest decrease at larger 

scales. The statistic is very stable; the six independent slices of the LCRS all yield consistent results. 
We apply the same statistical procedure to halos in numerical simulations of an open cosmological model 
and flat models with and without a cosmological constant. In contrast to the LCRS, all the models show 
a dispersion which increases for scales > \h~^ Mpc; it is uncertain whether this is a numerical artifact 
or a real physical effect. The standard cluster-normalized Cold Dark Matter model with Q.m = 1 as well 
as a tilted variant with n = 0.8 yield dispersions substantially hotter than the LCRS value, while models 
with low matter density — 0.3) are broadly consistent with the LCRS data. Using a filtered cosmic 
energy equation, we measure i7„i « 0.2, with small-scale bias factors b = 1.0-1.5 for high-density models 
and 6 — 0.7-1.1 for low-density models. 

Subject headings: cosmology — dark matter — galaxies: clustering — large-scale structure of universe 



1. INTRODUCTION 

The small-scale thermal energy of the observed galaxy 
distribution is an important diagnostic for cosmological 
models. For the past decade the pair velocity dispersion 
(712 (r) (Davis & Peebles 1983) has been the usual mea- 
sure of this quantity (e.g., Bean et al. 1983; de Lapparent, 
Geller, & Huchra 1988; Hale-Sutton et al. 1989; Mo, Jing, 
& Borner 1993; Zurek et al. 1994; Fisher et al. 1994; 
Marzke et al. 1995; Brainerd et al. 1996; Somerville, 
Primack, & Nolthenius 1997; Landy, Szalay, & Broad- 
hurst 1998; Jing, Mo, & Borner 1998). But in spite of its 
widespread application and the relative ease of its mea- 
surement within large redshift surveys, the ai2{r) statis- 
tic has a number of well-known deficiencies. Chief among 
them is its pair-wise weighting, which gives extreme in- 
fluence to rare, rich clusters of galaxies containing many 
close pairs with high velocity dispersion. 

Alternative statistics to measure the thermal energy dis- 
tribution have been suggested by Kepner, Summers, & 
Strauss (1997) and by Davis, Miller, & White (1997, here- 
after DMW). The Kepner et al. algorithm computes the 
pair-weighted dispersion as a function of the local galaxy 
density; this statistic demonstrates the heterogeneity of 
the environments of the local galaxy distribution, but it 
must be computed in volume-limited samples. The ci 
statistic described by DMW can be estimated within a 
flux-limited catalog and is readily interpreted in terms of a 



filtered version of the cosmic energy equation. The statis- 
tic is a measure of the rms one-dimensional velocity of 
galaxies, with large-scale bulk flow motions filtered out. 
DMW applied this statistic to the UGC catalog of opti- 
cal galaxies within the Optical Redshift Survey (Santiago 
et al. 1995), as well as the 1.2-Jy IRAS catalog (Fisher 
et al. 1995). They showed that = 1 simulations were 
far too hot to match the observed dispersion. Even when 
compared with simulations in which the small-scale kinetic 
energy had been artificially lowered by a factor of four, the 
observed velocity distribution was colder than the simu- 
lated distribution. 

However, the UGC catalog surveys a rather limited vol- 
ume of the local Universe, and the IRAS catalog is quite 
dilute and under-samples dense cluster regions. It is there- 
fore of considerable interest to apply the DMW statistic to 
a larger, more representative redshift survey such as the 
Las Campanas Redshift Survey (LCRS; Shectman et al. 
1996), and to compare the results with iV-body simula- 
tions of cosmological models which are favored by current 
data. This paper reports the application of this new statis- 
tic to the LCRS and compares the result to a few simu- 
lations of flat and open cosmological models. In a future 
paper (Baker, Davis, & Ferreira 1999), we discuss a wider 
variety of models, and we discuss in more detail the com- 
parison of the LCRS with Y-body simulations and the 
potential applications of cti as a cosmological probe. 
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2. APPLICATION OF cri TO THE LCRS 

The LCRS survey consists of 26,000 galaxies selected 
in a hybrid R band. The survey was conducted in six 
thin slices, each of size 1?5 x 80° on the sky, with median 
redshift cz = 30,000 km s'^ The redshift accuracy of 
the observations is typically (Terr = 67 km s~^ (Shectman 
et al. 1996), which is sufScient for measuring the thermal, 
small-scale velocity dispersion. 

For measurement of ui, we work with the subset of 
19,306 LCRS galaxies in the range 10,000 < cz < 
45,000 km s~^, and absolute magnitude —22.5 < M < 
— 18.5. To estimate the random background of the neigh- 
bors about eac;li galaxy, we used a catalog of 268,000 ran- 
domly distributed points with the same selection function 
as the LCRS galaxies, including the restriction against 
pairs with angular separation less than 55" caused by lim- 
itations of optical fiber placement. Since the six slices of 
the LCRS are spatially separated by more than the pro- 
jected separation used in the u\ statistic, the statistical 
procedure is applied to each slice individually and the re- 
sults are averaged. 

2.1. Method 

We now briefly describe our procedure, similar to that 
of DMW, for determining a\ . For each galaxy i in a slice 
of the survey, we lay down a cylinder centered on the 
galaxy in redshift space. Let r-p be the projected radius 
of the cylinder and vi its half-length along the redshift di- 
rection. For neighboring galaxies j within the cylinder, 
we construct the distribution Pi{Av), which counts the 
number of neighbors with redshift separation in a redshift 
bin centered at Av = Vj — Vi. The counts accumulated 
in Pi{Av) are weighted by the inverse selection function 
(pi/(j)j (though equal weighting yields virtually identical re- 
sults). We subtract from this distribution the background 
distribution Bi{Av), which counts the number of weighted 
neighbors expected for an unclustered galaxy distribution. 
We are interested in the width of the overall distribution 
D{Av) constructed by an appropriately weighted sum over 
the Ng galaxies: 



We achieve a more inclusive statistic by considering sep- 
arately the distributions of high- and low-density objects; 
that is, only galaxies with iVex > 1 are included in the sum 
for £)hi, while only galaxies with A'ex < 1 are included in 
the sum for Diq. Wc then weight the galaxies in the com- 
bined distribution according to 



-{ 



A^ex,i > 1 



< 1. 



(3) 



Here Aio and Ahi are normalization constants for the two 

distributions, chosen so that the distributions are weighted 
in proportion to the number of objects included: 



Nhi/Ng 



Ea. [Dhi{Av) - Dhi{oo)y 



(4) 



and similarly for ^iq. Here A^hi and Aio are the num- 
ber of galaxies with A'cx > 1 and A"cx < 1- respectively, 
thus A^hi + A^io = Ng. The baselines D{oo) are estimated 
from the flat tails of the distributions within 500 km s~^ of 
Av = ztvi. With this normalization the final distribution 
obeys D{Av) = 1. Note that scaling Dhi and Ao 
by the constants A does not affect the derived widths for 
these distributions; rather, it merely alters the weighting 
of the two in the combined distribution. 

This procedure, in contrast to that of DMW, allows 
us to include all of the available data, yielding an unbi- 
ased, object- weighted measure of the thermal energy of the 
galaxy distribution. It is the object- weighting which differ- 
entiates our procedure from the more traditional measure 
of the pair dispersion <Ti2{r); all galaxies (not pairs) are 
assigned equal weight in our statistic cti . 

We measure the width of the distribution D{Av) us- 
ing the convolution procedure outlined by DMW (equa- 
tion 18), in which a velocity broadening function f{v) is 
convolved with the two-point correlation function ^(r) to 
produce a model M{Av) = * / for D{Av): 



M{Av) 



j dr2-Kr j dy r"^ + y"^) f {Av 



y)- (5) 



1 



(1) 



where the weight for galaxy i is denoted by Wi. 

In order to make the statistic object- weighted rather 
than pair-weighted, we wish to normalize the distributions 
by the number of neighbors N^^ in excess of the random 
background, that is: 



W^ = Aex.i 



j2mM-B^iAv)] 



(2) 



This however presents a problem for galaxies which do not 
have enough neighbors to ensure that the sum is positive. 
DMW dealt with this problem by deleting these objects 

from consideration, but under half of the LCRS galaxies 
have at least one excess neighbor for Vp = lh~^ Mpc, and 
these galaxies are a biased sample because they populate 
over-dense regions. It is therefore desirable to modify the 
statistic to include galaxies with fewer neighbors. 



The two-point correlation function of the LCRS is well- 
approximated by ^(r) — (ro/r)'^, with ro = 5ft Mpc 
and 7 = 1.8 (Jing et al. 1998), while for the A^-body sim- 
ulations we use the cylindrically averaged mass correlation 
function (Av) measured directly from the particle dis- 
tribution. We find that an exponential broadening func- 
tion (see Diaferio & Geller 1996; Sheth 1996, Juszkiewicz, 
Fisher, & Szapudi 1998) 



/(.) = lexpf-M 



(6) 



provides a much better fit to the LCRS data and all A- 
body models than does a Gaussian. Here we have defined 
the width ai so that it is a measure of the rms velocity 
of individual galaxies in one dimension (with bulk motions 
on scales ;> lh~^ Mpc filtered out). The (object-weighted) 
rms difference in velocity between any two galaxies is then 
ai\/2 (DMW call this quantity, which is equal to the rms 
dispersion of the distribution /, the "intrinsic" dispersion 
CT/; we will work exclusively with cti to avoid confusion). 



Small-Scale Velocity Dispersion of the LCRS 



3 



0.10 



0.08 



0.06 



0.04 



0.02 



0.00 




-2000 -1000 1000 2000 

Av (km s"') 

Fig 1. — Galaxy-weighted velocity distribution D{Av) for 
the six LCRS slices. 

The three-dimensional dispersions are larger by an addi- 
tional factor -\/3- 

We perform a nonlinear ^^-minimization fit to deter- 
mine the width cri and amplitude of the model M(Av). 
Before fitting, we convolve the model with a Gaussian of 
rms (TcrrV^ = 95 km s~^ to account for the LCRS red- 
shift measurement uncertainties; the factor of ^2 converts 
from the measurement uncertainty for individual redshifts 
to the uncertainty for redshift differences, which are ac- 
cumulated in D{Av). We also include baseline terms in 
the model which are constant and linear in Aw, for a total 
of four fit parameters. The linear term is necessary for 
the LCRS because for simplicity we define "cylinders" in 
redshift space based on projected angular separation on 
the sky. This leads to a small gradient in the measured 
distribution function D{Av) because the "cylinders" are 
in fact conic sections, but the term is quite small because 
the length of the cylinders, 2vi, is much smaller than the 
typical redshift of galaxies in the survey. Although the 
gradient term has a negligible effect on the derived width, 
it does improve significantly the quality of the fit. 

2.2. Results for the LCRS 

We have used the six independent slices of the LCRS 
to estimate the errors in D(Av) as a function of Av in 
computing x^- However, we expect that the bins may be 
correlated due to sample variance; the fitting procedure is 
therefore not strictly legitimate, but the consistency of the 
results for the widths of the individual slices serves as a 
check on the degree to which sample variance affects the re- 
sult. We also expect > 1 if the exponential broadening 
function of width ai (assumed independent of r) provides 
an inadequate description of the small-scale velocities. 

The D{Av) distributions for the six independent LCRS 
slices are plotted in Figure ||, and Table [l| lists the de- 
rived widths. The second to last line gives the mean and 
standard deviation of the mean for separate fits to the 
six slices, while the last line is the result of a single fit 
to the combined distribution of all galaxies. Note that 



the dispersion measured for objects with excess neighbors 
{Ncx > 1) is clearly higher than that measured for objects 
with fewer neighbors. This behavior is expected because 
objects with more neighbors are found in regions of higher 
density, which tend to be hotter. 

The fit to the LCRS D{Av), shown in Figure ^, is quite 
good, with xl = 117/96 = 1.22; the probability of 
exceeding this value is 1 — P(x^|j^) = 7%. The best- 
fitting Gaussian f{v) is much worse, with xt = 1-84 and 



1 - P{x' 



10- 



Based on the mean of the six slices we adopt (Ti = 
126 ± 10 km s~^. This value has been computed for 
Tp = lh~^ Mpc and vi = 2500 km s~^. The results are 
quite insensitive to cylinder length, ranging only from 
117 ± 14 km s-i at vi = 1500 km s^^ to 132 ± 13 km s~i 
at vi = 3500 km s~^. Our chosen value vi = 2500 km 
is large enough to allow a clean measure of the tails of the 
distribution without significant non-linearities in the base- 
line gradient due to variations in the selection function. 

A modest decrease in cri is evident as rp is increased 
above Vp = Ih'^ Mpc (see Table |). AUhough the D{Av) 
distributions are very insensitive to rp, the averaged cor- 
relation function (Av) becomes broader as rp increases. 
As a result, smaller values of rp provide a cleaner measure 
of the true (real-space) velocity broadening on small scales, 
but decreasing rp below lh~^ Mpc reduces the signal-to- 
noise, as most galaxies have too few neighbors. The back- 
ground subtraction also becomes cleaner as rp is reduced. 

Note that for the larger value rp = 2h~^ Mpc used 
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Fig 2. — Velocity distribution D{Av) and fit for the com- 
bined LCRS data (upper panel), and residuals for the fit 
with errors estimated from the standard deviation of the 
six slices (lower panel; note the change in vertical scale). 
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Table 1 

The (Ti statistic for the six LCRS shces. The second to last line gives the mean and standard deviation of the mean of 
the independent slices, while the last line is the result of a fit to the entire dataset. We list cri for the combined dataset 
and for the subsets of galaxies with N^.,^ > 1 and A'ox < 1- N^-JNg is the fraction of galaxies with N^^ > 1, and Ng is the 
total number of galaxies. 



by DMW, our resuh is cti = 114 ± 10 km s^^. If, as 
in the DMW analysis, we do not account for broadening 
due to redshift measurement errors, the result increases to 
fJi — 136 ± 9 km s~^. Since the two surveys have compa- 
rable redshift uncertainties, our LCRS result is perfectly 
consistent with the value cri — 130 ± 15 km s^^ which 
DMW derived for the much smaller UGC catalog. 

3. COMPARISON TO AT-BODY MODELS 

We have completed a suite of A^-body simulations de- 
signed to predict the small-scale velocity dispersion in a 
variety of cosmological models. Here we discuss the results 
of a few of these models: the "standard" Cold Dark Mat- 
ter (SCDM) model and a tilted variant (TCDM), a model 
with a cosmological constant A (LCDM), and an open 
model (OCDM). The cosmological parameters for these 
models are listed in Table ^. All models are approximately 
normalized to the present-day abundance of clusters; the 
LCDM and TCDM models additionally satisfy the COBE 
normalization. The SCDM model is known to fail a num- 
ber of cosmological tests and is included for historical rea- 
sons, and only LCDM is fully consistent with current lim- 
its from high-redshift supernovae (Perlmutter et al. 1999). 
We note that on the scales relevant for our simulations, the 
TCDM power spectrum is indistinguishable from a rCDM 
spectrum with shape parameter F — 0.2. A broader range 
of models and a more detailed discussion of the simulations 
may be found elsewhere (Baker et al. 1999). 

Initial power spectra were obtained using the CMB- 
FAST code (Seljak & Zaldarriaga 1996, 1998). The sim- 
ulations were evolved on a 128^ mesh using a P M code 



(Brieu, Summers, & Ostriker 1995) in which short-range 
forces are computed using a special purpose GRAPE-3AF 
board (Okumura et al. 1993). We chose a box of size L — 
50h~^ Mpc to match the length of the LCRS cylinders; 
with Np = 64"^ particles this gives a mass resolution of 
1.3 X 10"rJ™/i-i Mq, where h = Hq/IOO km s'^ Mpc~\ 
A Plummer force softening e = 50h~^ kpc was used. The 
simulations were started at redshifts Zi = 15 (for = 1) 
or Zi = 19 (for il„i = 0.3) and evolved to z = in 1500 
time-steps using p = as the integration variable. 

The simulations are converted to "redshift" space by 
adding the velocities Vi along one of the three coordinates 
i to the positions xf. — > Xi -\-Vi/ H, where H is the Hub- 
ble constant. Periodic boundary conditions are applied at 
the box edges. We then apply exactly the same statistical 
procedure for determining ui as for the LCRS, except that 
the selection function is now unity. 

3.1. Tests of ai Measurements 

We have used our A^-body simulations to perform a 
number of checks on the robustness of our method for de- 
termining the small-scale velocity dispersion. One test is 
to ask how well our model is able to account for the red- 
shift measurement uncertainties in the LCRS. To simulate 
these uncertainties, we added Gaussian random velocities 
of rms (Terr along the "redshift" coordinate in the simu- 
lations. We then make two determinations of cti, which 
should ideally be equal. In one determination, the ran- 
dom velocities have been added and we perform an extra 
Gaussian convolution in the model to account for them. In 
the other, no random velocities are added and no Gaus- 



Tpjh-^ Mpc) cri(kms~^) A^hi/A^g 
0.5 136 ±10 0.23 

1 126 ±10 0.42 
1.5 107 ±8 0.55 

2 96 ±12 0.63 
2.5 99 ±13 0.68 

Table 2 

LCRS dispersion ai as a function of limiting projected radius Tp, and fraction of galaxies with excess neighbors. 
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Table 3 

Cosmological parameters for iV-body models: matter density fim, cosmological constant Ha, tilt n where P(k) cx fc", 
Hubble constant Hq — ft,/100 km s^^ Mpc~^, and rms mass fluctuation as in spheres of radius 8h~^ Mpc. 



sian convolution is necessary. We find that the two widths 
agree quite well, to within 10 km s""'^ over the range of in- 
terest for (Ti (100-300 km s~^). The agreement improves 
as (Ti increases and the uncertainties contribute relatively 
less to the width of the observed velocity distribution. 

A second test of the method is to compare velocity 
widths measured in real space with those measured in 
cylinders in redshift space. For this test, we replace the 
velocities of the simulation particles with velocities drawn 
from a random exponential distribution of a given rms a. 
It is straightforward to show that the velocity distribution 
appropriate for the difference distribution D{Av) is then 

Using this form in the redshift-space model (Equation S), 
we find that our procedure recovers the true velocity clis- 
persion with an accuracy better than 10% for ai in the 
range 100-300 km s^^. 



Finally, we can test the extent to which our measure- 
ment of CTi in the long redshift-space cylinders is contam- 
inated by motions on scales larger than lh~^ Mpc. First 
we construct distributions analogous to D(Av), but mea- 
sured in real space, with neighbors drawn from spheres of 
radius lh~^ Mpc in the simulations. These are compared 
to distributions with neighbors drawn from the long cylin- 
ders, also measured in real space. The widths of these 
distributions agree to within 1%, and we conclude that 
the contamination from large scales is negligible. 

3.2. Selection of Galaxies from the Mass Distribution 

We can easily compute ai for particles in the simula- 
tions, but the observed small-scale dispersion of galaxies, 
which correspond in some way to halos in the simulations, 
will in general differ from that of the mass. The internal 
velocity dispersions of galaxies are not included in the ob- 
served statistic; moreover, the galaxy population may be a 
biased tracer. In order to test whether our simulations can 
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Fig. 3. — Two-point correlation functions multiplied by (r/ro)^'^ for LCDM simulation halos for a range of Ns and a. 
Here tq = 5.1h~^ Mpc and the LCRS ^(r) appears as a solid horizontal line. In both plots, the curve at bottom (dotted 
line) shows Ng = oo (no halo splitting), and the solid curve shows the mass correlation function. In the plot at left, Ng 
is held fixed at 80, while a takes on the values 0, 0.25, 0.5 (top to bottom). In the plot at right, a is held fixed at 0.25, 
while Ns takes on the values 10, 20, 40, 80 (top to bottom). 
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Fig. 4. — Two-point correlation functions for the mass (dotted lines), halos (dashed), and LCRS (solid), plotted as in 
Fig. ||. Halos were selected using the parameters listed in Table 0. 



reproduce the LCRS result for (Ti, it is therefore impor- 
tant to identify "galaxies" within the A^-body simulations. 
Unfortunately the process of galaxy formation includes 
baryonic physics on a wide range of scales not probed by 
our dark-matter only simulations. For the present work, 
we define galaxies using a simple phenomenological model 
which we expect to yield similar results to those of larger 
gas-dynamical simulations. 

We first apply the standard friends-of-friends (FOF; 
Davis et al. 1985) algorithm to the simulations, with a 
linking length of 0.2 mesh cells and a minimum group 
size A^ > 10, corresponding to halos with mass M ^ 

lO^^rJm/i"^ Mq. We have also considered the HOP 
method (Eisenstein & Hut 1998) for defining halos, but 
we obtain similar results for reasonable parameter choices 
and do not discuss them here. 

Our limited resolution and the nature of the FOF al- 



gorithm lead to a serious and well-known over-merging 
problem, in which a large cluster containing many galax- 
ies will be identified as a single halo. This drastically 
lowers the small-scale velocity dispersion because the mo- 
tions of galaxies within clusters are neglected. To remedy 
this situation, we split halos with more than A^s particles 
by randomly selecting particles from within the halos and 
identifying these particles as galaxies. Halos identified in 
this way will again include the internal motions of galax- 
ies, but as the splitting is only applied to large, hot halos 
{Ns ^ 10), we expect these internal motions to have a neg- 
ligible effect on our result. Small halos with fewer than Ng 
particles are taken to be individual galaxies. 

For comparison with the LCRS, we wish to choose a set 
of halos which resemble the LCRS galaxies as closely as 
possible. Some A^-body models yield a correlation function 
^(r) which is too steep, and it is therefore advantageous 
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80 


0.00 
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310 
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40 
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80 
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Table 4 

Velocity width cti for the A^-body simulations. Results are listed for all particles and for the A^haios halos identified using 
our splitting procedure with parameters Ng and a. 



to select halos which are anti-biased on small scales (Jing 
et al. 1998). We accomplish this through our halo-splitting 
procedure by drawing random particles with a probability 
p which has a power-law dependence on the number of 
particles A'^ in the parent halo: p{N) — N"~^N~°', with 
a > 0. The number of galaxies per unit halo mass then 
falls as 7V~" for large halos. We choose parameters Ng and 
a which simultaneously mimic the power-law shape of the 
LCRS correlation function and produce approximately the 
correct number density of galaxies, n « 0.02h^ Mpc~^, im- 
plying 2500 galaxies per simulation volume. Increasing a 
tends to flatten the correlation function on small scales 
and yields fewer halos; increasing Ng at fixed a tends to 
lower the correlation amplitude and also yields fewer ha- 
los. This behavior is illustrated in Figure g|for the LCDM 
model. 

Figure ^ shows the correlation functions for our selected 
halos in each of the models. In the low-density models, we 
are able to select halos which match the LCRS ^(r) quite 
well. The normalization of the high-density models is such 
that ^(r) always falls below the LCRS power law on large 
scales. The TCDM halos match well at r <; 2h^^ Mpc. 
In the SCDM model, we are unable to reproduce exactly 
the shape of ^(r) without falling too far below the LCRS 
amplitude and prod ucin g too few halos. However, the dif- 
ferences in J2 (see §^^) computed from these correlation 
functions show that this mismatch should affect our esti- 
mate of fljn by at most 30%. 

3.3. Results for <ti 

The results for cti for our four cosmological models are 
listed in Table ^. We see that the mass in the two Vtm = 1 
models is far too hot on \h~^ Mpc scales, with cti well 
over twice the LCRS value. The spectral tilt of the TCDM 
model has very little effect on the small-scale velocities, as 
the result is nearly identical to the SCDM result. The 
mass in the low-fim models, on the other hand, is also 
hotter than the LCRS, but only by a factor of about 1.5. 

The halos in the simulations are somewhat cooler than 
the mass, with small-scale dispersions lower by factors 
in the range 0.7-0.9. The LCDM halos come closest to 
the LCRS value; at 143 km s^^, they are only marginally 
(1.7(t) hotter than the LCRS. The open model produces 
velocity dispersions slightly higher than the LCDM model, 
while the halos in the — 1 models are again much hot- 
ter than the LCRS data. 

Figure ^ shows that the exponential f(v) provides an 
excellent fit to the velocity distributions measured in the 
simulations in redshift space. We show distributions for 
the A^-body mass particles and for the halos. The halo 



distributions are noisier because there are many fewer ha- 
los than mass particles in the simulation volumes. The 
distributions for the SCDM and OCDM models are nearly 
indistinguishable from the TCDM and LCDM distribu- 
tions, respectively, and are not shown. 

We have also computed ui for galaxies drawn using more 
sophisticated semi-analytic techniques from a large Virgo 
simulation (Benson et al. 1999) of the LCDM model. This 
simulation has a mass resolution better than ours by about 
a factor of two, and the box length is nearly three times as 
large. The result is 126 km s^^, only slightly lower than 
our value of 143 km s~^. This suggests that our proce- 
dure for defining galaxies is reasonable. The Virgo re- 
sult exactly matches the LCRS dispersion, which suggests 
that the small-scale velocity dispersion predicted by the 
= 0.3 flat model is in fact perfectly consistent with 
the observational data. Further details of this comparison 
will be presente d in a future work (Baker et al. 1999). 

As noted in §2.2, the LCRS velocity width decreases 
somewhat as the limiting radius rp^,„ax is increased. In 
Figure |^, we show this scale dependence measured in in- 
dependent cylindrical shells of width \h~^ Mpc, where 
the limits on the radial integration in the model (Equa- 
tion ^) have been adjusted appropriately. Although the 
measured LCRS D{Av) shows little scale dependence, the 
integrated correlation function broadens with scale, lead- 
ing to a smaller measured velocity width. 

None of the A^-body models, however, are able to re- 
produce the scale dependence observed in the LCRS. The 
halos drawn from the Virgo simulation, which show very 
little scale dependence, come closest, while the other mod- 
els tend to show an increase in velocity dispersion with 
scale. Only the LCDM model is shown in Figure I but we 
find similar discrepancies for the other models as well. Al- 
though the rim = 0.3 LCDM model produces a reasonable 
match to the velocity dispersion on very small scales, all of 
the models seem unable to reproduce the observed coldness 
of the velocities on intermediate scales ~ l-3/i~^ Mpc. At 
present it is unclear whether this discrepancy is due to 
problems with the galaxy selection procedure, the resolu- 
tion of the simulations, or a more fundamental flaw in the 
cosmological models. 



3.4. Filtered Cosmic Energy Equation 

The (Ji statistic is ideally suited for the application of the 
cosmic energy (Layzer-Irvine) equation filtered on small 
scales. As shown by DMW, we expect a\ oc ^mJ2,m in 
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Fig. 5.— Velocity distributions (solid) and fits (dashed) for two representative A''-body models. In each plot the fit is 
shown as a dashed line and is generally invisible because it is obscured by the measured distribution. Top and bottom 
panels show the distributions for all particles and for halos, respectively. 



the absence of velocity bias, where 

J2= / ""''dr<(r). (8) 

The subscript m means that J2 is computed from ^m('"), 
the correlation function for the underlying mass. We can 
write this in terms of the measured ^(r) of an observed 
sample j by defining an effective bias 6^ = Ji,jl J2,m- If 
we then compare ctij measured for sample j with cti^jv 
measured for the underlying mass in an TV-body simula- 
tion with mass density parameter ri^r, we can measure the 
parameter 



If in addition we can choose a sample of A^-body halos 
which matches the correlation function of the sample j, 
then we have a direct measure of '^m'- 

f^rn = (crij/cri,Ar)^rijV, (10) 

where cti^at is now measured for the A^-body halos rather 
than the underlying mass. 

The results of combining the LCRS dispersion a\ = 
126 ± 10 km s~^ with our four cosmological A^-body mod- 
els are listed in Table ^. Based on the halos in each of the 
four simulations, we derive consistent values ~ 0.2. 
Note that the errors listed on are l-cr uncertainties de- 
rived solely from the LCRS a\ result; they do not include 
any systematic errors in the model results. The fact that 
we derive similar values of Sim from each of the different 
models is an important consistency check, and gives us 
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Model n^n ^m/b'^ b 

SCDM 0.22±0.03 0.10-0.14 1.2-1.5 

TCDM 0.18±0.03 0.12-0.15 1.0-1.3 

LCDM 0.23±0.04 0.20-0.27 0.8-1.1 

OCDM 0.19±0.03 0.20-0.29 0.7-1.0 

Table 5 

Density parameter and small-scale bias derived from the cosmic energy equation and the LCRS dispersion. 



confidence that our method is indeed a sensitive probe of 
the matter density. 

Table || also lists the values of flm/b'^ derived by com- 
paring the LCRS dispersion with the dispersion of the iV- 
body mass. The integral J2 converges rather slowly, and 
its value is quite sensitive to the integration limits rmin 
and Tniax- A reasonable lower limit is r„iin = 0.1h~^ Mpc, 
which eliminates from the analysis the internal velocity 
dispersion of typical galaxies and includes only the dis- 
persion of galaxies moving relative to each other. We 
might also take Tmax to be shghtly larger than lh~^ Mpc, 
since the length of the redshift-space cylinders means that 
there will be some contribution to cti from larger scales 
(although we have measured this effect in the simulations 
and have found that it is very small). The ranges shown for 
f2m/6^ were obtained by allowing and rmax to vary 
over the ranges 0.05-0.2 and 1-5 Mpc, respectively. 
Our results for the high-density models are consistent with 
the value fi^/fe^ = 0.14 ± 0.05 found by DMW, who only 
considered an flm = 1 model. 

The parameter ^Im/b'^ is approximately equal to /3^, 
where /3 ~ /b is the parameter measured by large-scale 
flow analyses. We find /? ~ 0.3-0.4 for the two high-density 
models, and (3 ~ 0.45-0.55 for the two low-density models. 
These ranges are generally consistent with some large-scale 
flow determinations (e.g., Willick & Strauss 1998; Baker 
et al. 1998; Davis, Nusser, & Wilhck 1996) but not with 
the POTENT analyses, which prefer /3 ~ 1 (e.g., Sigad 
et al. 1998). Of course, the bias may in general depend on 
scale, in which case our small-scale result need not match 
the (3 values measured using flows on much larger scales. 

Finally, we can combine the values of 0^ and flm/b'^ to 
obtain an estimate of the bias of the galaxy distribution 
on small scales. Our high-density models require biases 
b = 1.0-1.5, while the low-density models are slightly anti- 
biased, b — 0.7-1.1. These ranges are consistent with the 
biases measured directly from the correlation functions of 
the simulations. 

3.5. Effects of Streaming Velocities 

Although our goal is to measure the particle distribution 
function from redshift-space information alone, we must do 
this by considering the relative motions of pairs of galax- 
ies, for which we expect mean streaming as well as thermal 
motions. As defined in Equation ^, our model does not 
account for a non-zero first moment of the velocity distri- 
bution of pairs of galaxies. However, the first moment will, 
in general, be non-negligible due to the mean tendency of 
galaxies to approach each other, and it will contaminate 
a measurement of the second moment. On small scales in 
virialized clusters, for example, the infall velocity approxi- 
mately cancels the Hubble expansion, and so its presence 
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Fig 6. — Object-weighted velocity dispersion measured in 
independent cylindrical shells of width \h^^ Mpc. The 
LCRS data are shown as filled circles with error bars. Also 
shown are LCDM mass (squares) and halos drawn from 
our simulations (crosses) and from the Virgo simulation 
(triangles). 

can affect our measurements on lh~^ Mpc scales by of or- 
der 100 km s^^. Jing & Boerner (1998) have shown that 
the effect of the streaming motions on the estimate of the 
pairwise velocity dispersion can be dramatic, increasing 
ai2 from ~ 400 km s~^ to 580 km s^^ at \h^^ Mpc sepa- 
ration. 

The effects of the streaming motions can be incorpo- 
rated into our analysis by writing the distribution function 
in Equation ^ as 

/(«)=-^exp(-^), (11) 

where vi is the mean object- weighted streaming velocity, 
which is a function of separation, and u'l is the second 
moment of the streaming-corrected velocity distribution. 
The form of vi is unknown but can be measured directly 
from A^-body simulations. 

Our estimate of ci with iTf = will be smaller than a'l 
because streaming motions tend to cause objects to pile 
up at small velocity separations in redshift space. How- 
ever, (Ti has the advantage that it is a mo del- independent 
statistic, relying only on the assumption of an exponential 
velocity distribution. The comparison of the data with 
iV-body models is also consistent; to the extent that the 
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Table 6 

Velocity widths with corrections for self-similar streaming motions applied. Results are listed for the LCRS and for halos 
drawn from our LCDM simulation. 



models describe the real universe, the same streaming mo- 
tions will be present in both the data and the models, and 
will affect the estimates of <Ji similarly. Incorporating a 
non-zero vi introduces model dependencies into the mea- 
surement, and there is no guarantee that the infall mea- 
sured in the TV-body simulations matches that of the real 
universe. 

For the application of the cosmic energy equation, it 
is in fact more appropriate to use cti rather than a'l, be- 
cause contributions from both random thermal motions 
and mean streaming motions are already included. On 
the other hand, a'^ is a better measure of the truly ther- 
mal energy of the galaxy distribution. We can estimate it 
by using Equation |ll| with an appropriate model for vi. 
For the mean pairwise velocity, the simple form 



Vi2{r) = - 



FHor 



1 + (r/ro)2 - 



(12) 



(Davis & Peebles 1983) is often used, where F is a numer- 
ical factor, typically F = 1-1.5. Another expression has 
been proposed more recently by Juszkiewicz, Springel, & 
Durrer (1999): 



(13) 



where / 



n 



0.6 



1.2 — 0.6570 with 70 



-d\n^/dlnr\^^i, and 



[l+e(r)]e(r) = ^ / dxx'a^). 



(14) 



These two forms for V12 (r) are nearly equal at small scales 
r ^ lO/i-i Mpc if we set F = l-Sn^n^; note that F = 1 
corresponds to streaming motions which just cancel the 
Hubble expansion on small scales. 

Table ^ shows that the streaming correction has a sub- 
stantial effect on the derived LCRS velocity width, with a'l 
rising to 201 ± 13 km s^^ for = 1 and 261 ± 15 km s^^ 
for F = 1.8. The statistic worsens somewhat for 
F > 1. The A^-body models show similar behavior. We 
caution, however, that the streaming-corrected dispersions 
are model-dependent and are not an appropriate measure 
of the single-particle dispersion for use with the cosmic 
energy equation, which is defined in the comoving frame 
of the universe. This is in contrast to analyses of the pair 
dispersion, where it is appropriate to use the cosmic virial 
theorem, defined in the mean streaming frame. 

4. CONCLUSIONS 



Although the potential of small-scale cosmological ve- 
locities as a cosmological probe has long been recognized, 
the application of pair-weighted statistics is problematic. 
We apply an extended version of the more stable galaxy- 
weighted statistic of DMW to the Las Campanas Red- 
shift Survey. We derive a one-dimensional rms veloc- 
ity for individual galaxies relative to their neighbors of 
(Ji = 126 ± 10 km on scales ^ lh~^ Mpc. 

Using this new statistic, we find that the observed veloc- 
ities remain quite cold relative to the predictions of high- 
Qm A-body simulations. Tilting the power spectrum to re- 
duce the initial power on small scales does little to resolve 
this discrepancy. We have also examined flat and open 
models with fl„i = 0.3; these models produce significantly 
lower dispersions than the high-density models. Combin- 
ing the LCRS data with the predictions based on halos 
in the simulations, we measure consistent values fim ~ 0.2 
for all models, and we can rule out Jim = 1 with a high de- 
gree of confidence. Our result suggests that the extremely 
cold dispersion measured in the vicinity of the Local Group 
(Schlegel, Davis, & Summers 1994; Governato et al. 1997) 
might be a local anomaly, as currently popular low-density 
models can reproduce the observed mean dispersion on 
lh~^ Mpc scales. On the other hand, at slightly larger 
separations, we find evidence that all of the models may 
again be too hot relative to the observations. 

In the future, it will be extremely useful to apply our 
statistic to upcoming redshift surveys, such as the Sloan 
Digital Sky and 2dF surveys, which will contain enough 
galaxies to compute cri precisely for different sub-samples 
of the galaxy population. The Deep Extragalactic Probe 
(DEEP; Davis & Faber 1998) and other surveys at high 
redshift will also provide a measure of the evolution of 
(Ti, which can be used to place additional constraints on 
cosmological parameters and the bias of the galaxy distri- 
bution. 
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